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Abstract 

Using the fact that extremum of variation of generalized action can lead to the frac- 
tional dynamics in the case of systems with long-range interaction and long-term mem- 
ory function, we consider two different applications of the action principle: generalized 
Noether's theorem and Hamiltonian type equations. In the first case, we derive con- 
servation laws in the form of continuity equations that consist of fractional time-space 
derivatives. Among applications of these results, we consider a chain of coupled oscilla- 
tors with a power-wise memory function and power-wise interaction between oscillators. 
In the second case, we consider an example of fractional differential action 1-form and 
find the corresponding Hamiltonian type equations from the closed condition of the form. 
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1 Introduction 



Different physical phenomena such as anomalous transport or random walk with infinite mo- 
ments [U 12], dynamics of porous media [3], 0], continuous time random walk [3 El [7], chaotic 
dynamics [8] (see also reviews [9j [10]) can be described by equations with fractional integro- 
differentiation. Despite of fairly deep and comprehensive results in fractional calculus (see 
[TT1 fT2| [TBI [T4] ) a possibility of their applications to physics needs to develope specific phys- 
ical tools such as extension of fractional calculus to the areas as multi-dimension [Til [T7] . 
multi-scaling [151 HB], variational principles p^8l [19] . It became recently formalized the links of 
the dynamical equations with fractional derivatives to the systems with long-range interaction 
[20] |2TI [22] and long-term memory [23] . 

In this paper, we concentrate on two problems important for numerous physical applica- 
tions: conservation laws and Hamiltonian type equations, both obtained from the corresponding 
fractional action principles. In Sec. 2, we derive the Noether's theorem for a Lagrangian that 
includes non-local space-time densities. The Noether's theorem was also discussed in [261 EZ] • 
Our new derivation shows in an explicit way how fractional derivative in time emerges from 
the specific type of the memory function, and how fractional derivative in space is related to 
a specific long-distance potential of interaction (Sec. 3.) In Sec. 4, these results are applied 
to a chain of nonlinear oscillators that is a subject of great interest in statistics and dynamics 
[281129]. Finally, at Sec. 5, we derive a specific case of fractional Hamilton's equations. Different 
steps in this direction were performed in [321 1331 134"] . We consider the Lagrangian density as 
a functional without fractional derivatives but, instead, the differential 1-form has fractional 
differentials. Some examples are given for this type of systems. 

The main feature of this paper is the consideration of fractional type differentials or deriva- 
tives in both space-time coordinates. 
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2 Noether's theorem for long-range interaction and mem- 
ory 

2.1 Action and Lagrangian functionals 

Let us consider the action functional 



S[u] — d x d y jC(u(x) , u(y) , du(x) , du(y)) , (1) 
Jr Jr 

where x = (t,r), t is time, r is coordinate, and y = (t',r'), du(x) = (d t u(t,r),d r u(t,r)). The 
integration is carried out over a region R of the 2-dimensional space R 2 to which x belong. The 
field u(x) is defined in the region R of M. 2 . We assume that u(x) has partial derivatives 

dou{x) = du ^ r \ d lU {x) = dU ^ T , 

which are smooth functions with respect to time and coordinate. Here C(u(x),u(y), du(x), du(y)) 
is generalized density of Lagrangian. If 

C(u(x) , u{y) , du(x), du{y)) = C{u(x), du(x)) S(x — y), (2) 

then we have the usual action functional 

S[u] — / d 2 x C{u{x) , du{x)) . 
Jr 

The variation of the action ([I]) is 

5S[u,h}= [ d 2 x [ d 2 y (^—h(x) + — — rr d„h(x)+ 
Jr Jr \du{x) <9(<9 m m(x)) 

where fi = 0,1, = djdx^ and 

£ = C(u(x),u(y), du(x), du(y)), 
and h(x) = Su(x) is the variation of the field u. The variation ([3]) can be presented as 

5S[u > h] = J R £x i £y (Hy h{x) + w^T) d ^ x) ) ' (4) 



where 

C s = C(u(x),u(y), du(x), du(y)) + C(u(y), u(x), du(y), du(x)). 
We can define the functional 

1 r j2 „ 



L[x,u,du] = - d yC s , (5) 

1 JR 

which will be called the Lagrangian functional. For (jSJ), the functional (jSJ) is equal to the usual 
Lagrangian density C(u,du). 
Using notations 

SL[x, u, du] f l9 dC s 5L\x,u,du] f l9 dC s . . 

dy^r^i r/n / \ \ = / ^ 2/ "oTo 7 vT' ( 6 ) 



5u(x) J R du{xY 5{d^u(x)) J R d{d ll u{x)) 
the variation (jlj) is 

The structure of (J7J) is non-local, i.e., it includes, through (jSJ) and £ s a possibility of the long- 
term memory and long-range interaction. Let us derive the generalization of the Noether's 
theorem for this case in the way analogical to the local case (J2]). 

2.2 Equations of motion 

First, we separate the variation that is linked in the variation of coordinates 

x n _ x 'n = x n + ^ ( 8 ) 

and the variation caused by a change of the form of u, 

u(x) — > u'{x) = u(x) + 5u{x). 

The variation h(x) = 5u(x) of u(x) that is not due to the variation in coordinates is called 
local. The total variation is 

Au(x) = u'(x') - u(x) = 8u(x) + (du/dx^dx" (9) 
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in the first approximation with respect to Sx. Let us consider the variation of (CD) as 

5S[u,h] = [ d 2 x' L[x',u',d'u] - [ d 2 xL[x,u,du]. (10) 
Jr Jr 

The elements of the two-dimensional volume in new and old coordinates are related through 
the formula 

d 2 x' = J(x' /x)d 2 x, 

where J(x'/x) = det \dx' tl /dx u \ is the Jacobian of the transformation. Using the well-known 
relation 

det A = exp Tr In A, (11) 

and the linear approximation 

— = % + d v (6x*), 

we get 

J(x'/x) = l + d^Sx* 1 ). 

Here 8^ is the Kronecker symbol. 
For the variation ffTOl) . we get 

5S[u,h}= [ d 2 x {5L + L^(50) = 
Jr 

= [ d 2 x ( 5 -±6u + J±—5{d,u) + dpiLSx*)) , (12) 
where the variations of L are defined in (Q. Using 5(d^u) = d^(5u), and 

d,,5u = d„ ( — ^ ; Su I - d„ ( J L ; I 5m, 



we can rewrite (fl2l as 



+ / d 2 x^ (Mar" + J^Su J . (13) 
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The Gauss theorem gives 

f d 2 xd^ (l5x» + j77r-T$u) = [ dS^ (l6x» + J^5u) . (14) 

We assume that at the boundary of the domain of integration the function u(x) is selected 
in a definite manner (the boundary condition). Then stationary (a minimum or a saddle point) 
values of S[u] from variational equations 8S[u, h] — with 

[5u] dR = 0, [Sx] 8R = 

at the boundary OR of the domain of integration, constitutes the necessary and sufficient 
condition for the real evolution of the field, that is, that u = u(x) represents the true dynamics 
under the given boundary conditions. The stationary action principle gives 

5L[x,u,du] SL[x,u,du] _ 
5u(x) 6(d,u(x)) ~ 

where the variations of L are defined in ([6]). This equation is the Euler-Lagrange equation for 
Lagrangian functional L[x,u,du}. 

Let us consider three special cases of equation (}T5]) . 

(1) The absence of the memory and long-range interaction means that 

C(u(x),u(y), du(x), du(y)) = C(u(x), du(x)) 8(x — y). 

Then ( TT5|) gives the usual Euler-Lagrange equation 

dC(u,du) dC(u,du) _ 
du{x) "did^uix)) [ ' 

(2) If the generalized Lagrangian density is 

C(u(x), u(y), du(x), du(y)) = C(u(x), du(x)) c\{D, r) 8(x — y), 

where 

ci(D,r) = tL, (0<D<1), 

for a medium distributed on M 1 with the fractional Hausdorff dimension D, then Eq. (fl5l) has 
the form 

^ a ^-^^WM)-- a - (17) 
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This is Euler-Lagrange equation for the field u(x) = u(t, r) in fractal medium. Examples of the 
field (wave) equations for fractal medium string and fractional hydrodynamics are considered 
in [36]. For example, 

C(u{x), du{x)) = - (d t u{t, r)f - -v 2 (d r u{t, r)f , 

leads to the equation 

r) — v 2 d r (ci(D, r)d r u(t, r) \ = 

that describes the propagation waves in fractal medium. 
(3) Consider the action functional 

S[u] = J d 2 x J d 2 y (^d t u{x)g (x,y)d t ,u(y) - ^d r u(x)g 1 (x,y)d r >u(y) - V(u(x),u(y)) 

(18) 

where the kernels go(x,y) and gi(x,y) are responsible for nonlocal time-coordinates dynamics 
and V is nonlocal interaction potential. Then the Lagrangian functional (jSJ) is 

L[x, u, d t u, d r u] = 

d 2 y (\d t u(x)K (x,y)d t m(y) - \d T u(x)K 1 (x,y)d r iu(y) - U(u(x),u(y))) , (19) 



R 



where 



Ko{x,y) = ^[go{x,y) +g (y,x)), 

Ki(x,y) = ^\gi(x,y) +gi(y,x)], (20) 
U(u(x),u(y)) = V(u(x),u(y)) + V (u(y) , u(x)) . 



Assume that 



U(u(x),u(y)) = U(u(x))6(x - y). 
In this case, the Euler-Lagrange functional equation (|15p has the form (see also |23j ) 

f d 2 yd t K (x,y) d t ,u(y) ~ f d 2 yd r K x {x,y) d r ,u{y) + d ^M± = . (21) 
Jr Jr du{x) 
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It is an integro-differential equation, which allows us to derive field equations for different cases 
of the kernels K (x,y) and K\(x,y). 

In the absence of memory and for local interaction the kernels (1201) are defined at the only 
instant t and point r, i.e., 

K (x,y) = g 5(x - y), K x (x,y) = g x 5{x - y) 

with some constants go and g x . Then equation ( 12TI) gives 

g dMt,r)- gi dMt,r) + ^^=0. 

For example go = g x = 1, when 

U(u(t,r)) = — cos u(t, r), 

we get the sine-Gordon equation 

d 2 u{t, r) - d 2 u{t, r) + sin u(t, r) = 0. (22) 



2.3 Noether's current 



Let us derive the Noether's current by using the action variation ( fl3l) . 
The second integral of (TL3]) can be presented as 



d 2 xd u L5x^ + 



5L 



5{d lx u 



-8u 



d 2 x d„ 



5L 



-[5u+(d v u)5x u ] 



5L 



■(d v u)-6^L 



r M \5(d^u 

Using the total variation Eq. (|23l . and the energy-momentum tensor 

SL 



5x l 



6?, 



-{d v u)-5»L, 



the variation ffl3|) yields 



SS[u, h] 



d 2 x ( -J^- — d, 



5L . , 
n r,~ — - I du+ 
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(23) 



(24) 



5L 



+ I ^ x d » ( T7TT-J Au ~ 9 v 6x " ) ■ ( 25 ) 



Let us consider a continuous (topological) group of coordinate transformation x — > x' = 
x'(x, a), and let the field function u = u(x) admits the representation of this group: 

x — > x', u{x) — > u'(x'). (26) 

The invariance of functional S[u] with respect to ( |26l) means that 

/ d 2 x / L[x / ,u',du/dx']= / <i 2 a;L[a;, w, du/dx}. 
Jr Jr. 

The transformation ( 1261) constitutes a group. Therefore, infinitesimal forms of transformations 



026]) are 

Ax u = X^5a s , Au = Y s 5a s , {v = 0, 1, s = 1, .., m), (27) 

where and Y" s are the generators of the group of transformation correspondingly in the 
coordinate and the field representations. The index s — 1, .., m is defined by the representation 
of the group. 

For simple examples of transformation of the coordinate, time or scalar field, we have 

Ax" = X v 5a, (u = 0,l), (28) 

Au = Y5a. (29) 

Noether's theorem states that every continuous transformation of coordinate (128]) and field 
function (T291 . which ensures that the variation of the action is zero admits a conservation law 
in the form of a continuity equation [301 [31]. Substitution of ( 1281) . f l29l) . and (115!) into (123]) gives 

5%, /i] = /" d 2 x a M ( j^-tY - e^X u ) 5a. (30) 
Jr V^(^ n ) / 

In view of the fact that the variation of the parameter, 5a, is arbitrary, from (!30j) we get the 

conservation law 

= 0, (31) 

where 

J *=im y - KX "' (32) 
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is the Noether's currents and is defined in (I24j) . Equation ( 13TT) means that there exists 
conservation law. Nontriviality of equations fl3Tj) and fl32|) is that L and 6^ have nonlocal 
interaction and memory. 



3 Application of the Noether's theorem for long-range 
interaction and long-term memory 

3.1 Lagrangian functional and energy- momentum tensor 

For ii = 0, 1, x° — t, and x 1 = r, the Euler-Lagrange functional equation (ITST) is 

5L[x,u,d t u,d r u] d SL[x,u,d t u,d r u] d SL[x,u,du] 
Su(x) dt 5(d t u(x)) dr 5(d r u(x)) 

The time and space variables in action can be separated to consider the field with power-law 

memory and long-range interaction. Let K (x,y) and Ki(x,y) have the form 

K (x,y)=5(r-r , )JC (t,t'), (34) 

K 1 (x,y) = 5(t-t')JC 1 (r,r'), (35) 

and 

U(u(x),u(y)) = U(u(x))6(x-y), (36) 
where x = (t, r), and y = (t', r'). Then the Lagrangian functional is 

L[x,u,d t u,d r u] = - J dt'JCo(t, t')d t ud t >u — - J dr'K,i(r,r')d r ud r iu — U{u). (37) 

To present the long-term memory and long-range interaction, consider the kernels JCo(t,t') and 
/Ci(r, r') in the power-law forms 

£i(r,r') =Kx(\r-T , \) = — -, ^7 ri - t — T , (1 < a < 2), (38) 

U ' ; 1VI {J cos(7ra/2)r(2-a) (r-r'l"- 1 ' V ; ' V ; 



and 



M(t-t'), 0<t'<t; 
JCo(t,t) = { (39) 

0, a >t, f < 0, 



10 



with 



""-"- rd-ae-V (0</3<1) - (40) 



Then the variational derivatives of Lagrangian functional are 

6L dU(u(t,r)) 
5u du(t, r) 



(41) 



()L %D?u(t,r), (0 < /? < 1), (42) 



S(d t u 

5L 9 



•a— 1 



0/>w^=T u (*' r )' (K«<2), (43) 



5(d r u) d\r 
where f-Df is Caputo derivative [HI [12] defined by 

c<t^b , n 1 <^ r d n u(r,r) , _ n . 

9D?u(t,r) = —, -r -. ^ — -r — ^- L -^, (n- 1 < /3 < n). 

* v ' ; r(n - (3) J (t - Tf- n + l dr n v M ; 

Substitution of (j4"Tj) . fj42l and fT43l into (1331 gives the fractional field Euler-Lagrange equation 
^D? +1 u(t,r)-9^u(t,r) + = 0, (K a < 2, 0< < 1). (44) 



For a = 2 and (3 = 1, Eq. (1441) is the usual field equation 

#«(t, r) - ^(t, r) + = 0. (45) 

For the potential U(u(t,r)) = — cos u(t,r), equation (J4"5j) gives the sine-Gordon equation, 

<9 t 2 w(t, r) - d*u(t, r) + sin u(t, r) = 0, (46) 

and Eq. ( 1441) is a spatio-temporal fractional sine-Gordon equation 

Sf u - d^u + sin u = 0, (47) 

where we used for abreviation simplified notation for fractional derivatives. For the case (3 = 2, 
the equation was obtained in [20J. 

The energy-momentum tensor 6%, can be presented in the form 

9H = 95^ + t», (48) 
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where the first term is the diagonal part 

SL 



e = e° + e\ 



S(d t u] 



(dtu) + 



5L 



5(d r u 



■(d r u) -2L = U{u(t,r)). 



(49) 



This represents a pressure-like quantity. The second term in the rigt hand side of ( 1481) i 
nondiagonal of 9%: 



is a 



5L 



T . 



-(d v u) 



5L 



where 



r o _ 1 _ Q 
'o ~~ 'i ~~ u > 



J^(d r u) = [ ^ii(f,r)] d r u(t,r), 



5L 



r o = 



5(d r u 



-td t u) = -g 



— u(t,r) dtu(t,r). 



d\r 



(50) 

(51) 
(52) 

(53) 



For r e M 3 the spatial components of r/f, represent shear stress tensor. The value 9 represents 
a pressure-like quantity, normal stress. 



3.2 Homogeneity in time 

The homogeneity in time means invariance of action with respect to the transformation 

t — > t' = t + a, r — > r, u —> u. 
Then the infinitesimal transformations are 

Ax» = 5%5a, Au = 0, (/i = 0, 1) 

with the generators 

X» = 6», Y = 0, (/i = 0,1). 
The Noether's current has two components 



5L 



6(du,u) 



d t u - 5%L , (jm = 0, 1) 



(54) 



(55) 



(56) 



(57) 



Using x° = t, and x l = r, we get the continuity equation 



d t J° + drJ 1 = 0, 



(5f 
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where 



J° = - (j^j (d t u) -L)=- (d t u(t, r) f Df«(t, r) — L 



= ~ (\d t u(t,r) $D?u(t,r) + I^ r «(t,r)^p-«(t,r) +Uj , (59) 

and 

SL d a ~ l 
As a result, the continuity equation is 

d t (±d t u c ^ t u + 4^^=™ + U )~9r {ddtu-^u^ = 0, (61) 
where u = u(t,r). Here 

H = \d t u °D p t u + ±gd r u-^-^u + U (62) 
is a fractional generalization of density of energy (density of Hamiltonian) , and 

ga-l 

V = -9 d tUg^ziu (63) 

is a fractional generalization of density of momentum. For the case a = 2 and f3 = 1, we obtain 
the well-known relations [3T] : 

n = \{d t uf + \g{d r uf + U, V = -gd t ud r u. (64) 

The continuity equation ( IBTj) can be presented in an usual form of the conservation of energy 

d t H + d r V = (65) 

with the fractional generalizations of the energy and momentum given by ( 1621) and ( 1631) . 



3.3 Homogeneity of space 

The homogeneity of space means invariance of action with respect to transformations 

r — > r' = r + a, £ — > £, u —>■ u. (66) 
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The corresponding infinitesimal transformations are presented by 



Ax" = 5%5a, Au = (67) 

with the generators 

X" = 5%, Y = 0. (68) 
Then the Noether's current has the following two components 

= -01 = _ ( J^^u) - 5?L ) , = 0,1). (69) 



Using x° = t, and x 1 = r, we get the continuity equation that corresponds to the homo- 
geneity of one- dimensional space 

d t J° + d r J l = 0, (70) 

where 

J° = -J^( 9 r u ) = -d r u(t,r) %D?u(t,r), (71) 

and 

Jl = ~s^uj^ + L = gdrU ^ r) W\^ u ^ r) + L = 

= \dtu{t,r) $Dlu(t,r) + ±gd r u(t,r)^^u(t,r) - U. (72) 
As a result, the continuity equation (1701 is 

- dt (d r u °D?u) + d r (±d t u °Dlu + ^ g d r u-^^u - u) = (73) 

that can be interpreted as the momentum conservation law in the case of fractional dynamics. 
For (3 = 1 and a = 2, equation (173!) gives 

- d t (d r ud t u) + d r Q(3tu) 2 + \g{d r uf - Uj = 0. (74) 

Note that, in general, for integer and fractional values of a and (3, J® 7^ Jq, and the energy- 
momentum tensor 6% is not a symmetric with respect to \x and v. 
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3.4 Field translation invar iance 



For the case U(u) = 0, the action (|18p is invariant with respect to transformations 

x fj> x 'fi = a; M ) it — > it + a. 

The infinitesimal transformations are presented as 

Aa^ = 0, Au = 5a 

with the generators 

X^ = 0, F = l. 

The Noether's current 

J" = T^— V (A* = 1,2) 

has components 

J° = ^y= ?I2fu(t,r), (0 </?<!), 
As a result, the continuity equation 

fyj* 4 = d t J° + drJ 1 = 

has the form 

d t $D?u(t,r) - gd r ^—- x u{t,r) = 0, (1< a < 2, < (5 < 1). 
For a = 2 and /3 = 1, we get 

dtu{t,r) -gd 2 r u{t,r) = 0, 

which is the usual field equation ( Pf5i) for [/ (w) = 0. 
Note that [H] 

3 f23f«(t,r) = %Dl +(} u(t,r) + * _ d t u(0, r). 
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To have the relation d t QD^u(t,r) = qD\ ; tt(t, r), the initial conditions d t u(0,r) = should 
be applied. In general, equation flHTj) cannot be presented as 



^D^u(t, r) - 9jjrrzu(t, r) = 0, (1< a < 2, < (3 < 1), 



34) 



which is the fractional field equation (I44p for C/(w) = 0. 

This shows that the conservation law flHTj) . in general, doesn't coincide with the field equation 
as it is happened for integer derivatives, unless, we use special boundary and initial conditions. 

4 Chain with long-range interaction 



It was shown in [201 Ell [221 [2H [25] how the long-range interaction between different oscillators 
can be described by the fractional differential equations in the continuous medium limit. In 
this section, the Noether's theorem will be applied to such kind of systems. 



4.1 Equation of motion and Noether's currents 

Let us define the action as 



S[u n ] 



4£ 



Kit) - KM*)) 



+oo 

U{u n {t),u m {t))), (85) 



n,m=— oo 



where u n is displacement of the n-th oscillator from the equilibrium, 

U(u n (t)) = -g Q J a (\n - m\) (u n (t) - u m (t)) 2 



and 



J a (\n-m\) = - -— , (a>0). 

\n — m\ a+1 



The Lagrangian of the chain is 



i ul(t) - V(u,(t)) 



j +oo 

7^0 ^a(l^ - "^D (M„(t) - M m (t)) 2 



n,m=— oo 



(86) 



*7) 



(88) 
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The equation of motion 

_d£__l^£_ = Q 

du n (t) dtdii n (t) v 1 

for Lagrangian (1551) has the form 



dp OV ( ~] 

— u n (t)+ ™ +g Ja{\n-m\)[u m {t)-u n {t)\ = 0. (90) 



dt 2 riW duJt) 

v ' n,m=— oo 



A continuous limit of equation (1901) can be defined by a transform operation from w n (t) to 
u(x,t) [20, EH E21 EU EE]. First, define w n (t) as Fourier coefficients of some function u(k,t), 
k e [-K/2,K/2], i.e. 

+oo 

u(t,k) = u n {t)e- ikx -= T A {u n {t)}, (91) 

n=— oo 

where x n = nAx, and Ax = 2ix / K is a distance between nearest particles in the chain, and 

i r+K/2 

u n (t) = — dk u(t, k) e ikXn = Fl l {u{t, k)}. (92) 

K J-K/2 

Secondly, in the limit Ax — > (K — > oo) replace u n (t) = (2tt/ K)u(x n , t) — > u(x,t)dx, and 
x n = nAx = 27m/ K — » x. In this limit, Eqs. (19~iT) . (I9"2"j) are transformed into the integrals 

/+oo 
rfx e- ikx u{t,x) = F{u{t,x)} = lim T A {u n {t)}, (93) 
-oo Aas '^° 

i p+oo 

u(t, x) = — dk e lkx u(t, k) = f- l {u(t, k)} = lim F7}{u{t, k)}. (94) 
Applying ( |9TI) to ( |90l) and performing the limit ( l93l) . we obtain 

where 

/ ixa\ 

Qot = 2g Q {Ax) a T{-a) cos [-) (96) 

is the renormalized constant. Equations (I9l)j) were considered in [2"Ul I2"lj I2"2"l |2"H [23]. 
Consider a continuous transformation of time and field 

t - w n (t) - <(*'), (97) 
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which form a continuous group with generators X and Y n such that the infinitesimal transfor- 
mations are 

At = X5a, Au n = YJa. (98) 
The corresponding Noether's current is 

+°° f)r 



where 



-oo 



8f 



is the energy. Then 

+00 

du n (t) 



dL r 1 
J= 2 + £X ( 10 °) 



The equation of conservation law is 

d 

— J = 0. 101 
In the next subsections, we consider examples of the conservation laws. 



4.2 Homogeneity of time 

The homogeneity of time means the invariance of action with respect to the transformation 
(compare to ( 1541) ) 

t — > t + a, u n — > u n . (102) 

Its infinitesimal form is 

At = 5a, Au n = 0, (103) 

with generators 

X = l, Y n = 0. (104) 

Then the Noether's current is 

+ OO Q£ +OO 

j t = -e = - E wm <(Jt)+£ = ~ 52 < &® +£ = 
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+oo ^ +00 ^ +00 

E 9 £ ^W*))"Tffo ^ J a (\n-m\)(u n (t)-u m (t)) 2 = -H, 



7i=— co n=— 00 n,m=— 00 



where if is the Hamiltonian. The conservation law df/y dt — is in the continuous limit gives 
the equation 



d r 00 , 

dr 



X - [dtu(t, r)f + V (u(t, r)) + ^0, w(t, r) 



d 
dt 



+00 

drH = 0, (105) 



dt 

where 1 < a < 2, and the density TC of Hamiltonian fl62|) is introduced. 

One can compare this equation with equation ( 1651) for (3 = 1. Both results coincide if the 
boundary conditions 



lim P = (106) 

r— >±oo 



is applied. 



4.3 Translation invariance 

For V(u n ) = 0, the action (185]) is invariant with respect to the transformations 

t — > t, u n — > u n + a, (107) 



or in the infinitesimal form 



with the generators 



Then the Noether's current is 



At = 0, Au n = 8a, (108) 



X = 0, Y n = 1. (109) 



where P is the tolal momentum. The conservation law is dP/dt = 0. 
The continuous limit of this conservation law gives the equation 



d_ 

dt. 



drd t u(t,r) = 0. (Ill) 
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Let us compare this equation with Eq. ( IHTl) that is derived for scalar field u(t, r). Integration 
of (|8ip with respect to coordinate r gives 



/+oo /-+oo <9 Q 

dr°D?u(t,r)-g y drd r -^ ZT u(t,r) = 0, (1< a < 2, < /? < 1). (112) 

Then 

/ + °° / (R a— 1 \ +oo 

dr °D?u(t,r) - g{^—u(t,r)) ^ = 0. (113) 

As a result, equation flHTT) coincides with fillip , for /3 = 1, if we use the boundary conditions 

lim J?— U (t,r) = 0. (114) 

r— >±oo O r 

5 Fractional Hamilton's equations 
5.1 De Donder-Weyl Hamiltonian 

The idea of variation of fields based on a manifestly covariant version of the Hamiltonian 
formalism in field theory known in the calculus of variation of multiple integrals [391 BOj has 
been proposed by M. Born and H. Weyl [H] . The mathematical study of geometrical structures 
underlying the related aspects of the calculus of variations and classical field theory has been 
undertaken recently by many authors (see for example H21 H3]). 

The De Donder-Weyl Hamiltonian form of the field equations [39] are 

where 

dC 

d(d^u) 

is called the multi-momenta, and 

H(u, tt") = (dpu) ^ - C (117) 

is called the De Donder-Weyl Hamiltonian density function, C = C(u,dfj,u) is Lagrangian 
density. These equations are known to be equivalent to the Euler-Lagrange field equations if C 
is regular in the sense that 

r d 2 C 1 

^ 0. (118) 



det 



d{d lx u)d(d tl u) 
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5.2 Hamilton's equations of integer order 

Let us consider Hamiltonian systems in the extended phase space of coordinates u, 71"^). 
The evolution of fields is defined by stationary states of the action functional 

S[ u , tt] = J [tt^u - H(u, ?r)] d 2 x, (119) 

where the Hamiltonian density 7i is defined by (I117p . both u and it are assumed to be indepen- 
dent functions of x = (t, r). In classical field theory, the evolution of the field u(x) is derived by 
finding the condition for which the action integral (11191) is stationary (a minimum or a saddle 
point). The action functional (11191) can be rewritten as 



S[u,tt} = j (120) 

where 

cu fl = 7i> 1 du-n(u,7i)dx> 1 . (121) 

is the Poincare-Cartan 1-form or the action 1-form. 
The stationary action condition SS = leads to 

= 0. (122) 

Here the exterior derivative is 

d = {dx v )D x u + (du)D u + (d7r")iV, (123) 
where we put new notation for the derivatives 

D " = ^ n " = l' c "" = a|- < 124 > 

The De Donder-Weyl Hamiltonian form of field equations in fl 122[) can be obtained from 
the condition that the derivative is exterior one. This condition is equivalent to the stationary 
action principle 5S[u, n] = 0. Then Eq. ( 1121ft gives 

duo^ = d(ir^du) - d(Hdx^) = (D x v-n^)dx v Adu+ (D u Ti^)du A du + (D 7T ,n >1 )d-K u A du- 
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- (D x »H)dx u A dx" - {D u H)du A dx^ - {D^H)d-K v A dx 11 . (125) 
Using dn Adx = —dx A dn, and (D x ^7i), D u n = 0, we get 

duo = (D x tt + D u 7i) dx A du — (D n ndu — D n 7idx) A dir, (126) 

where the matrix form of the equation is used. The relation D^n = 1 gives 

duj = (D x 7t + D U H) dx A du — (rfw — D n Hdx) A c?7r. (127) 

From f)122p . we have 

du - D w Hdx = 0, D x ii = -D U H. (128) 

As the result, we obtain 

°" D^H, ^ = -D u H, (129) 



dx^ dx^ 
which are the well-known De Donder-Weyl Hamilton's equations ( 11151) . 

5.3 Fractional Hamilton's equations 

The fractional generalization of the form (11211) can be defined [23 ES] by 

u(a) = irdy-H(u,ir)d°x. (130) 

It will be called the fractional Poincare-Cartan 1-form or simply the fractional action 1-form. 
We can consider the fractional exterior derivative of the form (11301) . and use 

d a u(a) = (131) 

to obtain the fractional field equations. Here the fractional exterior derivative is 

d a = d a s x v D« v +dyD° + dy u D«„, d a s x = T" 1 (a + l)d a x a , (132) 



where D% v , D", D" v can be fractional derivatives of different types [Hj. For example, for 
x^ G 1R 2 , such that x = (t, r), we use 



D2 = (°D?>,B?) t 
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where f.D" is Caputo fractional derivative, and D® 1 is the Riesz derivative. Fractional differ- 



ential forms have been suggested in [37] and it is used to describe dynamical systems [35], [38] 
Then, by some transformations (see Appendix), one can obtain 



d a uj(a) = [L>> + DIH\ d a 8 x A d a s u - 



71 



-rdZu-DZHdZx 
1 (2 — a) 



Using ffT33D and flTBTj) . we get 

.l-a 



71 



, . -dy - Dpid%x = 0, ££tt = —D"H. 
r(2 - a) 



For the case u = u(x) 



d a a u = d a s xDy. 



As the result, we obtain 



(133) 



(134) 



(135) 



L>> = (7^)^(2 " <*)D^H, = -D a u U. 



(136) 



These equations are the fractional generalization of De Donder-Weyl Hamilton's equations 

dm. 

As an example, we can consider 



n(u, tO = VV) 2 - ^A- 1 ) 2 + u(u). 



(137) 



Using 



D ^ f _ r(g + i: 



r(/? + i-a 

and r(z + 1) = zT(z), equations fl 136[) give 



(138) 



Df°U = 7T°, D^U = -7T 1 , 



(139) 



A a °7ro + A Ql vri = -^f/H 
After substitution of (11391) into (11401) . we obtain 



(140) 



(D?°) 2 u(t,r) - (D^) 2 u{t,r) + D^U(u) = 0. 



(141) 
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For qj = ot\ — a u — 1, Eq. ( 11411) is the usual wave equation 

d 2 u(t, r) - d 2 .u(t, r) + 



du(t, r) 

For U (u) = 0, Eq. ( 11411) has the form 

(D?°) 2 u(t,r) - (D^) 2 u(t,r) = 0. (142) 

The solution u(t, r) of this equation is a linear combination of the solutions U\(t, r) and u^it, r) 
of the equations 

D?° Ul (t, r) - D^ Ul (t, r) = 0, D?°u 2 (t, r) + D^u 2 (t, r) = 0. (143) 

For «i = 1, there exists a relation between the Dirac solutions and the fractional extension of 
D'Alembert expression that is considered in |44j . 
Using the property 

D«D a x u = D 2 «u, (144) 

and we get the fractional equation 

D 2ao u - D 2ai u + D**U(u) = 0. (145) 



For a special case, 



equation (I145j) gives 



r(2) i+ Qti 



U(u)=m 2 y u 1+a % (146) 
1 (2 + a u ) 



D 2ao u{t, r) - D 2ai u{t, r) + m 2 u(t, r) = 0. (147) 

This is a fractional generalization of Klein-Gordon equation (see also [201 EH])- Note that 
Eq. (I147p is not Lorentz invariant equations. To obtain fractional relativistic equations, the 
fractional power of D'Alembertian should be used [15] (see also Sec. 28 of [H]). The causality 
principle [31], HS] also should be taken into account. 
For 

U(u) =sm(u(t,r)) (148) 
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equation (11451) gives 

D 2a °u - D 2 r ai u + sin(u + a u vr/2) = 0. (149) 

This equation is fractional sine-Gordon equation that was considered in [20] for cto — 1 and 
a„ = 2. 

6 Conclusion 

Exploiting a variational principle to obtain fractional dynamics seems to be a fairly powerful 
tool that permits a universal consideration of situations with fractal time and space. In this 
paper, we have demonstrate it by deriving different equations of using generalized Noether's 
theorem. Such equations are similar to the regular conservation laws, although the presence 
of fractional time derivatives reflects a dissipation. Similar variation of action can be used to 
derive Hamiltonian type equations although the situation here is not uniquely defined and some 
freedom of choosing the differential action 1-form leaves different possibilities. 

Acknowledge 

This work was supported by the Office of Naval Research, Grant No. N000 14-02- 1-0056, and 
the NSF Grant No. DMS-0417800. 

References 

[1] E.W. Montroll, M.F. Shlesinger, "The wonderful world of random walks" In: Studies in 
Statistical Mechanics, Vol. 11. J. Lebowitz, E. Montroll (Eds.), (North-Holland, Amster- 
dam, 1984) pp.1-121. 

[2] J. P. Bouchaud, A. Georges, "Anomalous diffusion in disordered media-statistical mecha- 
nisms, models and physical applications" Phys. Rep. 195 (1990) 127-293. 



25 



[3] R.R. Nigmatullin, "The generalized fractals and statistical properties of the pore-space of 
the sedimentary-rocks" Phys. Stat. Sol. B 153 (1989) 49-57; "Fractional integral and its 
physical interpretation" Theor. Math. Phys. 90 (1992) 242-251. 

[4] R. Gorenflo, F. Mainardi, "Fractional calculus: integral and differential equations of frac- 
tional order", in Fractals and Fractional Calculus in Continuum Mechanics (Springer, 
Wien, 1997), pp.223-276. 

[5] V.V. Afanasiev, R.Z. Sagdeev, G.M. Zaslavsky, "Chaotic jets with multifractal space-time 
random walk " Chaos 1 (1991) 143-159. 

[6] R. Metzler, J. Klafter, "The random walk's guide to anomalous diffusion: a fractional 
dynamics approach" Phys. Rep. 339 (2000) 1-77. 

[7] E. Barkai, R. Metzler, J. Klafter, "From continuous time random walks to the fractional 
Fokker-Planck equation", Phys. Rev. E 61 (2000) 132-138. 

[8] G.M. Zaslavsky, "Chaos, fractional kinetics, and anomalous transport" Phys. Rep. 371 
(2002) 461-580. 

[9] R. Hilfer (Ed.), Applications of Fractional Calculus in Physics, (World Scientific, Singa- 
pore, 2000). 

[10] B. West, M. Bologna, P. Grigolini, Physics of Fractal Operators (Springer, New York, 
2003). 

[11] S.G. Samko, A. A. Kilbas, O.I. Marichev, Fractional Integrals and Derivatives Theory and 
Applications (Gordon and Breach, New York, 1993). 

[12] I. Podlubny, Fractional Differential Equations (Academic Press, San Diego, 1999). 

[13] K. Miller, B. Ross, An Introduction to the Fractional Calculus and Fractional Differential 
Equations (Wiley, New York, 1993). 

[14] A. A. Kilbas, H.M. Srivastava, J.J. Trujillo, Theory and Application of Fractional Differ- 
ential Equations (Elsevier, Amsterdam, 2006). 

[15] M.M. Meerschaert, D.A. Benson, B. Baeumer, "Operator Levy Motion and multiscaling 
anomalous diffusion", Phys. Rev. E 63 (2001) 1112-1117. 

[16] G.M. Zaslavsky, "Multifractional kinetics" Physica A 288 (2000) 431-443. 



26 



[17] M.M. Meerschaert, J. Mortensen, S.W. Wheatcraft, "Fractional vector calculus for frac- 
tional advectiondispersion" Physica A 367 (2006) 181-190. 

[18] O.P. Agrawal "Formulation of Euler-Lagrange equations for fractional variational prob- 
lems" J. Math. Anal. Appl. 272 (2002) 368-379; 

D. Baleanu, O.P. Agrawal, "Fractional Hamilton formalism within Caputo's derivative" 

Czechoslovak J. Phys. 56 (2006) 1087-1092. 
[19] R. Schumer, D.A. Benson, M.M. Meerschaert, S.W. Wheatcraft, "Eulerian derivation of 

the fractional advection-dispersion equation" J. Contaminant Hydrol. 48 (2001) 6988. 
[20] N. Laskin, CM. Zaslavsky, "Nonlinear fractional dynamics on a lattice with long-range 



interactions" Physica A 368 (2006) 38-54 pjnTS I/0512010^. 
[21] V.E. Tarasov, CM. Zaslavsky, "Fractional dynamics of coupled oscillators with long-range 
interaction" Chaos 16 (2006) 023110; "Fractional dynamics of systems with long-range 



interaction" Commun. Nonlin. Sci. Numer. Simul. 11 (2006) 885-898 ( |dEEP S/0512013). 
[22] N. Korabel, CM. Zaslavsky, V.E. Tarasov, "Coupled oscillators with power-law interaction 
and their fractional dynamics analogues" Commun. Nonlin. Sci. Numer. Simul. accepted 
( |math-ph/0603074p . 

[23] V.E. Tarasov, CM. Zaslavsky "Fractional dynamics of systems with long-range space 



interaction and temporal memory" Physica A (2007). Preprint math-ph/0702065 



[24] N. Korabel, CM. Zaslavsky, "Transition to chaos in discrete nonlinear Schrodinger equa- 



tion with long-range interaction" Preprint math-ph/0607030 



[25] V.E. Tarasov, "Continuous limit of discrete systems with long-range interaction" J. Phys. 

A 39 (2006) 14895-14910. 
[26] J. Cresson, "Fractional embedding of differential operators and Lagrangian systems" 



Preprint |math.D S/0605752 



[27] G.S.F. Frederico, D.F.M. Torres, "A formulation of Noethers theorem for fractional prob- 



lems of the calculus of variations" Preprint math.OC/0701187 



[28] J. Frohlich, R. Israel, E.H. Lieb, B. Simon, "Phase transitions and reflection positivity. I. 
General theory and long range lattice models " Commun Math Phys 62 (1978) 1-34; 



27 



H. Brezis, E.H. Lieb, "Long-range atomic potentials in Thomas- Fermi theory" Commun. 
Math. Phys. 65 (1979) 231-246. 
[29] J. Barre, F. Bouchet, T. Dauxois, S. Ruffo, "Large deviation techniques applied to systems 
with long-range interactions " J. Stat. Phys. 119 (2005) 677-713; 

J. Barre, D. Mukamel, S. Ruffo, "Inequivalence of ensembles in a system with long-range 
interactions" Phys. Rev. Lett. 87 (2001) 030601. 



[30] E. Noether, "Invariant variation problems" physics/0503066 M. A. Tavels English trans- 
lation of "Invariante Variationsprobleme" , Nachr. d. Konig. Gesellsch. d. Wiss. zu Got- 
tingen, Math-phys. Klasse, (1918) 235-257; Transport Theory and Statistical Physics 1 
(1971) 183-207. 

[31] N.N. Bogoliubov, D.V. Shirkov, Introduction to the Theory of Quantized Fields (Wiley, 
New York, 1980). 

[32] F. Riewe, "Nonconservative Lagrangian and Hamiltonian mechanics" Phys. Rev. E 53 
(1996) 1890-1899; "Mechanics with fractional derivatives" Phys. Rev. E 55 (1997) 3581- 
3592. 

[33] D. Baleanu, S.I. Muslih, "Lagrangian formulation of classical fields within Riemann- 



Liouville fractional derivatives" Phys. Scripta 72 (2006) 119-121 Qhep-th/05 10071] ); 
S.I. Muslih D. Baleanu, E. Rabei, "Hamiltonian formulation of classical fields within 
Riemann-Liouville fractional derivatives" Phys. Scripta 73 (2006) 436-438; 
E.M. Rabei, K.I. Nawafleh, R.S. Hijjawi, S.I. Muslih, D. Baleanu, "The Hamilton formal- 
ism with fractional derivatives" J. Math. Anal. Appl. 327 (2007) 891-897. 
[34] A. A. Stanislavsky, "Hamiltonian formalism of fractional systems" European Phys. J. B 49 
(2006) 93-101. 

[35] V.E. Tarasov "Fractional variations for dynamical systems: Hamilton and Lagrange ap- 



proaches" J. Phys. A 39 (2006) 8409-8425 dmath-ph/ 0606048). 
[36] V.E. Tarasov, "Continuous medium model for fractal media" Phys. Lett. A 336 (2005) 
167-174 ( |cond-m at / 0506 137[ ) ; " Fractional hydrodynamic equations for fractal media" Ann. 
Phys. 318 (2005) 286-307 dphysics/0602~096| ). 



28 



[37] K. Cottrill-Shepherd, M. Naber, "Fractional differential forms" J. Math. Phys. 



42 (2001) 2203-2212 (math-ph/0301013); "Fractional differential forms II" Preprint 



math-ph/0301016 



[38] V.E. Tarasov, "Fractional generalization of gradient and Hamiltonian systems" J. Phys. 

A 38 (2005) 5929-5943; "Fractional generalization of gradient systems" Lett. Math. Phys. 

73 (2005) 49-58; "Fractional statistical mechanics" Chaos 16 (2006) 033108. 
[39] Th. De Donder, Theorie invariantive du calcul des variations (Nuov. ed. Gauthier-Villars, 

Paris, 1935); 

H. Weyl, "Geodesic fields in the calculus of variations" Ann. Math. 36 (1935) 607-29. 
[40] H.A. Kastrup, "Canonical theories of Lagrangian dynamical systems in physics" Phys. 

Rep. 101 (1983) 1-167 and references therein. 
[41] M. Born, "On the quantum theory of the electromagnetic field" Proc. Roy. Soc. (London) 

A 143 (1934) 410-437; 

H. Weyl, "Observations on Hilbert's independence theorem and Born's quantization of 

field equations" Phys. Rev. 46 (1934) 505-508. 
[42] G. Giachetta, L. Mangiarotti, G. Sardanashvily, New Lagrangian and Hamiltonian Methods 

in Field Theory (World Scientific, Singapore 1997). 
[43] I.V. Kanatchikov, "Canonical structure of classical field theory in the polymomentum 



phase space" Rep. Math. Phys. 41 (1998) 49-90 ( hep-th/ 9709229). 

[44] T. Pierantozzi, L. Vazquez, "An interpolation between the wave and diffusion equations 
through the fractional evolution equations Dirac like" J. Math. Phys. 46 (2005) 113512. 

[45] C.G. Bollini, J.J. Giambiagi, "Arbitrary powers of DAlembertians and the Huygens prin- 
ciple" J. Math. Phys. 34 (1993) 610-621. 

[46] S. Kempfle, "Causality criteria for solutions of linear fractional differential equations" 
Fractional Calculus and Applied Analysis 1(4) (1998) 351-364; 

S. Kempfle, I. Schafer, H. Beyer, "Fractional calculus via functional calculus: Theory and 
applications" Nonlinear Dynamics 29 (2002) 99-127. 



29 



Appendix 

To prove the proposition (j!33p . we use the rule 



00 fi s 



and the relation [H] 
for integer s, where 
For example, we have 



dx s 

s=0 



<V [d>] = ( 8 >1) 



dx s 

_ {-lf- l aY{k - a) 
" k) r(l - a)T(k + 1) ' 



00 as 
s=0 ^ ' 

= d a s x u A t*X («) (££„A M ) = (D^A^) dX A rfX- 

As a result, 

rf Q c;(a) = d a (irdy) - d a (Hd a s x) = 

= (Dy)d°x a dy + {Dy)dy a g?> + (Dy)dy a d>- 

- {piH)d a s x A t£ac - {DlH)d a s u A d> - (D*H)dy A d>. (150) 

Here _D" = f Q D^ is Caputo fractional derivative, and = d a /d\r\ a is the Riesz fractional 
derivative [H]. Using 

dy A g?> = -d^x A d> 

and D x 7i(u, n) = 0, -D"7r = for Riesz and Caputo derivatives, we can rewrite equation (11501) 
in the form 

d a uj( a ) = [Dy + D a u n\ dy a dy - [(Dy)dy - (D%n)dy\ a dy. (151) 

Substitution of 

7T 1 -" 



D > = T{2^y (152) 



into equation (11511) gives (11331) . 
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